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Abstract 

The generalized Breit-Wheeler process, i.e. the emission of e~^e~ pairs off a probe photon prop- 
agating through a polarized short-pulsed electromagnetic (e.g. laser) wave field, is analyzed. We 
show that the production probability is determined by the interplay of two dynamical effects. The 
first one is related to the shape and duration of the pulse and the second one is the non-linear 
dynamics of the interaction of with the strong electromagnetic field. The first effect manifests 
itself most clearly in the weak-field regime, where the small field intensity is compensated by the 
rapid variation of the electromagnetic field in a limited space-time region, which intensifies the 
few-photon events and can enhance the production probability by orders of magnitude compared 
to an infinitely long pulse. Therefore, short pulses may be considered as a powerful amplifier. The 
non-linear dynamics in the multi-photon Breit-Wheeler regime plays a decisive role at large field 
intensities, where effects of the pulse shape and duration are less important. In the transition 
regime, both effects must be taken into account simultaneously. We provide suitable expressions 
for the e~^e~ production probability for kinematic regions which can be used in transport codes. 

PACS numbers: 13.35.Bv, 13.40.Ks, 14.60.Ef 
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I. INTRODUCTION 

The rapidly progressing laser technology developments ^ offer opportunities for inves- 
tigations of quantum systems with short and/or intense pulses Several fundamental 
processes of electron-photon interactions in the nonlinear regime thus become accessible. 
Once these are understood experimentally and theoretically, e.g. within the framework of 
the standard model of particle physics or plain quantum electrodynamics (QED), one can 
search for new phenomena hinting also to "new physics". Among the elementary electro- 
magnetic (e.m.) interaction processes is the "conversion of light to matter". Generically, 
this notion refers to the emergence of particles coupling to the e.m. field. Having in mind 
electrons (e~) and positrons (e"^) one is interested in the conversion rate into e^, their phase 
space distributions, the back-reaction on the original e.m. field etc. 

Several variants of such conversion processes are known. The linear Breit- Wheeler pro- 
cess jsl 7' + 7 — )■ e+ + e~ refers to a perturbative QED process; the generalization to the 
multi-photon process 7' -|- 717 — )■ e"*" -|- e~ (nonlinear Breit-Wheeler process) were done in 
the pioneering work of Reiss j4| as well as Narozhny, Nikishov and Ritus js], Attributing 
theses processes to colliding null fields one can imagine another aspect. In the anti-node 
of suitably counter propagating e.m. waves an oscillating purely electric field can give rise 
to the dynamical Schwinger effect j^j; in the low- frequency limit one recovers the famous 
Schwinger effect js] awaiting still its experimental verification. These kinds of pair creation 
processes are related to highly non-perturbative effects jol, [lo|. Once pair production is 
seeded in very intense fields further avalanche like particle production can set in which then 



could screen the original field or even limit the attainable field strength One can re- 
late the Breit-Wheeler process to the absorptive part of the probe-photon correlator in an 
external e.m. field; in our case the latter being a null field too. 

In the present paper we focus on colliding null fields in the multi-photon regime and 
consider the generalized Breit-Wheeler effect for short pulses of e.m. wave fields ranging 
from weak to high intensities. Phrased differently we analyze e^e~ pair production by a 
probe photon 7' traversing a coherent e.m. (i.e. laser) field. The latter one is characterized 
by the reduced strength 

e 

where {A^) is the mean square of the e.m. potential, and is the electron mass (we use 



natural units with c = h = 1, e^/An = a 
variable characterizing both null fields is 



1/137.036). A second relevant dimensionless 
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where Sthr = 4Mg and s = 2a;a;'(l — cos 6^p) (for head-on collision geometry, 9^^, = vr); u' 
and k, k' are the frequencies and three-wave vectors of the laser field and the probe photon, 
respectively. The variable C is a pure kinematic quantity with the meaning that for C > 1 the 
linear Breit- Wheeler process 7' + 7 — ?■ e"^ + e~ is sub-threshold, i.e. kinematically forbidden. 
However, multi-photon effects enable the non-linear process 7' + 727 — + e~ even for 
C > 1 which we refer to as sub-threshold pair production. The non-linear Breit- Wheeler 
process has been experimentally verified in the experiment E-144 at SLAG |l2|. There, the 
minimum number of photons involved in one e~^e~ event can be estimated by the integer 
part of C(l + 'C^), i-e., five. (To arrive at such an estimate recall that the reduced strength ^ 
is related to the laser intensity II via ~ 0.56(cj(eV))~^10~"'^^lL/(W/cm^), and therefore, 
at u' = 29 GeV, u = 2.35 eV, and at peak focused laser intensity of 1.3 x 10^®W/cm^, 
one gets ^ = 0.36 and ( = 3.83. The laser pulses contained about thousand cycles in a 
shot, allowing to neglect the details of the pulse shape and duration.) A laser intensity of 
~ 2 X 10^^ W/cm^ has been already achieved 13 [.Intensities of the order of II ~ 10^'^. ..10^^ 
W/cm^ are envisaged in near future at the CLF 1^, ELI 15[, and HiPER 16[ laser facilities. 
Such large laser intensities allow for larger values of ~ II compared to the SLAG E-144 
experiment. 

The new generations of optical laser beams are expected to be essentially realized in 
short pulses (with femtoseconds duration) with only a few oscillations of the e.m. field. 
High laser intensities are presently achieved by the chirped pulse amplification resulting 
in short pulses. As shown for the Gompton effect in 17H22I] and for the Breit- Wheeler 



effect in 
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261] the pulse shape and the pulse duration become important. That means the 



treatment of the intense laser field as an infinitely long wave train is no longer adequate. 
Keeping the spatial plane-wave character we are going to explore here the e~^e" production 
as generalized Breit- Wheeler process in finite pulse approximation (FPA), i.e. investigate the 
impact of the temporal pulse structure, and provide the conditions under which the infinitely 



long pulse approximation (IPA) can be applied. This problem is of practical 



investigation of e~^e production in transport Monte Garlo calculations 
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interest for the 



28| . where the 
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probability of pair production in a background field is taken as an input. 

We show below that the e^e~ production probability is determined by the non-trivial 
interplay of two dynamical effects. The first one is related to the shape and duration of the 
pulse. The second one is the non-linear dynamics of the in the strong electromagnetic 
field, independently of the pulse geometry. These two effects play quite different roles in 
two limiting cases: The pulse shape effects manifest most clearly in the weak-field regime 
characterized by small values of the product ^ (. The rapid variation of the e.m. field in very 
short (sub cycle) pulses enhances strongly few-photon events such that their probability may 
exceed the IPA prediction by orders of magnitude. Non-linear multi-photon dynamics of the 
strong electromagnetic field plays a dominant role at large values of In the transition 
region, i.e. at intermediate values C,"^ ~ 1, the pair creation probability is determined by the 
interplay of both effects which must be taken into account simultaneously. 

Our paper is organized as follows. In Sect. II we derive the basic expressions for the 
probability of e^e~ creation in FPA and consider a few prototypical pulse envelope shapes. In 
Sect. Ill we discuss the case of ultra-short (sub cycle) pulses where the number of oscillations 
of the laser field is smaller than one. The case of short pulses with a few oscillations of 
the laser field within one pulse is considered in Sect. IV. In particular, we analyze the 
enhancement of the production probability in the sub-threshold region at small values of 
discuss the case of intermediate ~ 1, and evaluate the production probability at large 
values of C,'^. Our conclusions are given in Sect. V. In Appendix A, for completeness and 
easy reference, we present some details of the derivation of the production probability for 
very high intensities, 1. 

II. ELECTRON-POSITRON EMISSION IN A SHORT PULSE 
A. General formalism 

In the following we employ the e.m. four-potential of the circularly polarized laser field 
in the axial gauge A'^ = (0, A{(f))) with 




(3) 



where (p = k ■ x is invariant phase with four- wave vector k = (w, k), obeying the null field 
property k^ = k ■ k = {a dot between four- vectors indicates the Lorentz scalar product); (p 
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is the carrier envelope phase; |a. 



7t |2 



a^, d^ciy = 0. Transversahty means ka^^y = in 



the present gauge. Instead switching on/off the periodic e.m. field we encode the finiteness 
of a pulse in the envelope function /(0) with lim /(0) = (FPA). To characterize the 
pulse duration one may use the number of cycles in a pulse, N = A/n = \tuj, where the 
dimensionless quantity A or the duration of the pulse r are further useful measures. Below 
we analyze the dependence of observables on the shape of /(</>) for a variety of relevant 
envelopes. The IPA case is defined by = 1. The carrier envelope phase is particularly 
important if it is comparable with the pulse duration A. In IPA it is anyhow irrelevant; in 
FPA with ^ A the production probability would be determined by an involved interplay 



of the carrier phase, the pu 



emphasized, e.g., in 
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se duration and pulse shape as well as the parameters ^ and C, as 



30|. In present work, we drop the carrier phase, thus assuming 



^ A, and concentrate on the dependence of the production probability on the parameters 
^ and C together with pulse shape and pulse duration. A detailed analysis of the impact of 
(p on the pair production needs a separate investigation which is postponed to subsequent 
work. 

Utilization of the e.m. potential of ([3]) leads to two significant modifications of^the tran- 
sition amplitude in FPA compared to IPA. In IPA, the Volkov solutions 
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32l | refer to 



Fermions with quasi-momenta and dressed masses. In FPA, all in- and out- momenta and 
masses take their vacuum values. The finite (in space-time) e.m. potential ([3]) for FPA 
requires the use of Fourier integrals for invariant amplitudes, instead of Fourier series which 
are employed in IPA. The partial harmonics become thus continuously in FPA. The S matrix 
element is expressed generically as 



S 



—le 



fi 



dl Mfi{l){2ny6\k' + lk-p-p') 



(4) 



where k, k', p and p' refer to the four-momenta of the background (laser) field ([3]), incoming 



probe photon, outgoing positron and electron, respectively^ 



similarly to the case of the non-linear Compton effect 
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he transition matrix Mfi{l), 



2l| . consists of four terms 



Mfi{l) = ^M«C«(/) , 



(5) 



i=0 
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where 



2vr J 

— oo 
oo 

C<(i)(/) = — / rf0/2(0)e*^^-*^(<^) , 
27r J 

— oo 
oo 

— oo 
oo 

C'^'\l) = ^ I d<Pf{<P) sin0e^'<^-^W , (6) 

— oo 



with 



V{(t>) = z j cos(0' - 0o)/(0') - e'C^ / d(^'f\(t>'). (7) 

— oo — oo 

The quantity z is related to ^, I, u = {k' ■ kY/ {4:{k ■ p){k ■ p')), and ui = l/( via 



z = 21^,1^(1-^] . (8) 

The phase 0o is equal to the azimuthal angle of the electron emission direction in the e^e~ 
pair rest frame (j)p> and is related to the azimuthal angle of the positron as (po = (f)p + vr. 
Similarly to IPA, it can be determined through invariants ai^2 as cos 0o = c^i/^, sin (pQ = a2/z 
with ai^2 = e (ai^2 ■ p/k ■ p — ai^2 ■ p'/k ■ p'). 

The transition operators M^*-* in Eq. ([5]) have the form 

M(^) = Up. M^^ Vp (9) 

with 

^ ' 4(A;-p)(A;-p') 

2(A;-p') 2(A;-p) ' 2(A: ■ p) 2(A; ■ p') ' 

where m and are Dirac spinors of the electron and positron, respectively, and e' is the 
polarization four-vector of the probe photon. 



6 



The integrand of the function C^^^ does not contain the envelope function fjd)) and 



therefore it is divergent. One can regularize it by using the prescription of Ref. 17|. The 
formal result 

oo 

C(°HO = ^ / d<f^{^^os{<P-<Po)f{<P)-eCufi<P)) e^'^-^^W + e-^^(°) (11) 

— oo 

contains a singular term (last term) which however does not contribute because of kinemat- 
ical considerations implying / > 0. The differential probability of e'^e~ pair production in 
terms of the transition matrix Mfi{l) in Eq. (jl]) reads 

oo 

c 

It may be represented in conventional form as a function of u and 0p 

dW aMeC/^ 1 



dcppdu 167riVo v?l'^\/u — 1 



dl wil) (13) 



C 



with 



\w{l) = (2«, + l)|C(°)(/)p + e'(2«-l)(|C(2)(0P + |C(3)(/)|2) 

+ ReC(°)(/) (e'C«(0 - ^(aiC^^HO + «2C(3H0))* • (14) 

The differential probability dW in Eq. f ll3p . in fact is the probability per unit time (or 
rate). The time units in IPA (AT^^^^)) and FPA (AT^^^^)) are different. The ratio of 
Nq = AT^^^^y AT^^^"^^ may be evaluated as following. The variation of e.m. energy of a 
pulse in a volume 5*^2;, where S is an unit cross section in the x — y plane, per AT^^^^^ 
is equal to the integral of the energy flux vector for the electromagnetic energy I (Poynting 
vector) over the area dS I. Taking into account that this integral is the same in FPA and 
IPA one finds 



/ Sdz{E%p^ + B 



FPA) 



JSdziEjp^ + Bjp^) 



-oo 





where A = 27r/a; is the wave length, and E = — ^ and B = V x A are electric and magnetic 
fields, respectively. For a convenient comparison of IPA and FPA results, the latter one is 
scaled in Eq. ^ by I/Nq. For IPA, No = l. 



B. Envelope functions 



We consider one-parameter and two-parameter envelope functions. Among the one- 



parameter functions we choose the hyperbohc secant (hs) and Gaussian (G) pulses 

/hs(0) = — ]-j , /g(0) = exp 
cosh 
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2ra, 



(16) 



As the two-parameter function, we choose the symmetrized Fermi (sF) shape 33 1 



/ (0) ^ + (17) 

cosh ^ -|- cosh I 

The scale parameters A, tq and Tsf determine the normalization factor Nq in ( !T5|) : A^hs = 
e (1 + si-) ' = ^ (l + 4) , and A^^p = ^ (f, (t) log ^^^^^ + (t)) with t = 
(cosh^ + l)/sinh^. In the latter case we have defined Fi{t) = {t^ + l){-t'^ + lOt^ - 
l)/16t and F2{t) = (3t^° - 35t® + 90t^ - 90t^ + 35t'^ - 3) /(24(t2 - 1)3). in the limit of small 
6/rsF ^ 0, one finds N,f = ^ (l - f;^) + O (exp [-f^]) . 

For large A and Tsf, and small b/Tsp ^ 1, one can find Agp — A^hs at Tsf = A. Therefore, 
for the sake of comparison we denote A^o = A^hs, Tsf = A and choose the ratio b/A as 
an independent parameter, which determines in turn the normalization factor Nsf at finite 
values of b/A. 

The scale parameters A and tq are related to each other by tq = y/nNo ^1 H 
for the fixed normalization factors Nq = Nq = A^^s- 



2ttN§ 



The one- and two-parameter envelope functions are exhibited at the left and right panels 
of Fig. [U Top panels depict to the envelope functions /(0), the middle and bottom panels 
depict the product /(0)cos0, which determines the function V in Eq. ([7j). The top and 
the middle panels are for an ultra-short pulse (sub cycle) with the number of oscillations 
less than one, A^ = 0.5. The bottom panels correspond to a short pulse with A^ = 5. One- 
parameter envelopes are similar to each other and are close to the two-parameter envelope 
with 6/A ^ 0.5. Decreasing b/A results in an essential modification of the envelope function 
/(0): it becomes close to the flat-top profile with the double step shape ^(A^ — 0^). 

III. ULTRA-SHORT PULSES 

In this section we consider the pair production due to interaction of the probe photon 
with an ultra-short pulse, where the number of cycles less than one. 
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FIG. 1: Pulse envelope f{^^ (top panels) and the product /(</>) cos(^ (middle and bottom panels) 
as a function of the invariant phase ^. The left and right panels correspond to the one- and two- 
parameter envelope functions, respectively. The top and the middle panels exhibit an ultra-short 
pulse with the number of oscillations less than one, = 0.5, while the bottom panels are for a 
short pulse with N = 5. 

A. The case of small field intensity (^^ ^ 1) 

Consider first the case of small field intensities and a finite sub-threshold parameter ( 
characterized by the relations z 1 oi ^( <^ 1. 

The basic functions C(^)(/) in Eqs. (E]) and can be expressed in this regime as a 
superposition of the functions 



oo oo 

3^(0 = ^ / d<f>^'"^fi<f>)9{<P)= I dqF{l-q)G{q) 



(18) 



where F{p) and G{q) are the Fourier transforms of the envelope function /(0) and the func- 

oo oo 

tion^7(</.) = exp [-zP(</.)], respectively: F{p) = ^ / #e*^"^ /(</.), G{q) = ^ J d<pe<'^ g{<P). 

— oo — oo 

For small values of 2;, 2; <^ 1, G{q) ^ 5{q — go); where go — ("P^) with go ~ ^ 1; ^-iid 
~ F{1). Keeping the leading terms in Eq. ( fT4l) with C^^^ ~ 3^(/ — 1) ~ F(/ — 1), one 
can obtain an approximate expression for the total production probability: 

00 

W = aMeC^^^^ J dl^l) F\l - 1) , (19) 
C 

with 

1 ^ 

$(/) = v ! dcosO (--'^ + u-^ , (20) 
J \ui uf 2J 



where u = 1/(1 — t;^cos^6'); 6 and v are the polar angle and the velocity of the outgoing 



positron in the e^e c.m.s., respectively: v = a/1 — C/^- An explicit calculation results in 

The Fourier transforms of the envelope functions ( !T6|) and (fT7|) read 

A 



^"^^^^ " 2coshi7rA/ 

^g(0 = ^=exp 
V 27r 



2 



F..(0- ; + ^-Tt; . (22) 

1 - exp [-f ] smhTrM 

The square of the Fourier transforms of the envelope functions for a sub cycle pulse with 
= 0.5 are presented in Fig. O On the left panel, the solid and dashed curves correspond 
to the hyperbolic secant and Gaussian shapes, respectively. One can see a fast monotonic 
decrease of both the functions with some enhancement of at large values of /. The square 
of the Fourier transform for the symmetrized Fermi shape is shown in the right panel, where 
the solid, dot-dashed and dashed curves correspond to the ratios b/A = 0.1, 0.3, and 0.5, 
respectively. One can see large qualitative and quantitative differences between the one- 
parameter and symmetrized Fermi shapes, in particularly, at 6/A < 0.3. In the second case, 
decreases exponentially as exp [— 27rA-^]. The slope decreases proportionally to b/A (at 
fixed A). Also, the function oscillates with the period 61 = vr/A = n/O.bn = 2. Contrary 
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FIG. 2: Square of the Fourier transforms of envelope functions for a sub cycle pulse with N = 
0.5. Left panel: The solid and dashed curves depict the hyperbolic secant and Gaussian shapes, 
respectively. Right panel: The solid, dot-dashed and dashed curves show the symmetrized Fermi 
shape for b/A = 0.1, 0.3, and 0.5, respectively. 

to the above one-parameter shapes, the function Fsf has a significant high-/ component at 
2 < / < 4. This strong effect is not seen in the 0-space (cf. Fig. [H top panels), where all 
envelope functions look similar to each other. However, the differences in /-space are very 
important for the pair production. 
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FIG. 3: The probability W of e^e~ production as a function of the sub-threshold parameter 
^ for one-parameter envelope functions for an ultra-short pulse with N = 0.5. The solid and 
dashed curves correspond to numerical calculations with the hyperbolic secant and Gaussian shapes, 
respectively. The symbols "star" and "plus" are for the approximation (|19p . The thin solid 
curves marked by dots correspond to IPA. The left and right panels are for = 10"^ and 10~^, 
respectively. 



Our prediction for the total probability of e"'"e~ pair production as a function of the 
sub-threshold parameter ( for the one-parameter envelope functions for an ultra-short pulse 
with = 0.5 is shown in Fig. [31 The solid and dashed curves exhibit results of numerical 
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calculations using Eq. (fT3|) with the hyperbolic secant and Gaussian shapes, respectively. 
The symbols "star" and "plus" display the resultsobtained by using the approximation (fT9|) . 
The thin solid curves marked by dots correspond to the IPA case. The left and right panels 
display results for ^"^ = 10^^ and 10~^, respectively. One can see an agreement of predictions 
for the ultra-short pulse and IPA near and above the threshold at C ^ 1) ^^^d a strong 
difference between them below the threshold, i.e. for C > 1- Our approximate (analytical) 
solution of Eq. f|T9l) is in a fairly well agreement with the complete numerical calculation. 
The function $(/) in Eq. (fT9|) is rather smooth compared to the Fourier transform F{1 — 1), 
therefore, the dominant contribution to the integral in Eq. flTI?]) comes from the lower limit 
of / and, qualitatively, the slope of the probability as a function of C, is determined by the 
scale parameters of the envelope functions 

WUO ~ exp [-ttAC] , Wg{0 ~ exp [-rl^] . (23) 

Despite of the exponential decrease of the probability W as & function of one can see a 
large difference (several orders of magnitude) between predictions for the ultra-short pulse 
and IPA (or "crossed field approximation"). In the latter case the probability decreases 
much faster with increasing Q. 




FIG. 4: The same as in Fig. [3] but for symmetrized Fermi shape envelope. The solid, dot-dashed 
and dashed are for 6/A = 0.1, 0.3 and 0.5, respectively. The corresponding approximate solutions 
are shown by symbols "+", "x" and respectively. 

Our results for the symmetrized Fermi envelope is presented in Fig. |H Now, the shape 
of the probability is determined by the two parameters h (or 6/A) and A 

sin^ AC . (24) 



W^sf(C) ~ exp 



-2.A-C 
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The first term describes the slope of the probability as a function of (. The slope is propor- 
tional to the "ramping time" of the envelope function, b (or to the ratio h/ Delta at fixed 
A). The second term, following from the Fourier transform shown in Fig. HI describes some 
oscillations with a period inversely proportional to the duration A of the fiat-top section; 
it is independent of the ramping parameter h. Again, one can see a great difference be- 
tween predictions for the ultra-short pulse and IPA on qualitative and quantitative levels. 
The probability in IPA has a typical step-like behavior, where each new step indicates the 
contribution of the next integer harmonic. In FPA, the probability decreases monotonically 
with a slope determined by the shape of the envelope. The quantitative difference is rather 
large and, as predicted by results shown in Figs. [3] and HI can reach orders of magnitude 
depending on the shape of the envelope(s). 

B. Anisotropy 

As we have shown above, at small values of ^, 2; <^ 1, the probability of e^e~ production 
is essentially determined by the pulse shape. The function g{(j)) in Eq. ( [T8|) is not important 
and, therefore, the total probability would be isotropic with respect to the azimuthal angle 
0e- = 00 because only the function 7^(0) contains a 0o-dependence. For finite values of 2, the 
function g{(j)) becomes important, and the electron (positron) azimuthal angle distribution 
is anisotropic relative to the direction of the vector in Eq. 02]), at least for the monotoni- 
cally rapidly decreasing one-parameter envelope shapes. The reason of an anisotropy is the 
following. At finite values of z, the function Y{1) in Eq. f lTSjl is determined by the integral 
over dcj) with a rapidly oscillating function proportional to the exponential 



In case of a fast-decreasing function f{(f)'), the contribution of the term proportional to 
sin 00 is much smaller compared to the term proportional to cos 0o. At finite z, the dominant 
contribution to the functions y{l) comes from the region where the difference in the exponent 
is minimal, i.e. 0e = 0o — 0. This means that the electrons would be emitted mostly along 
the vector a^^ and the positrons in the opposite direction. 
We define the anisotropy of the electron emission by 




(25) 



A 



dW{(P,) - dW{(Pe + tt) 
dW{(j)e) + dW{(j)e + 7r) ■ 



(26) 
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FIG. 5: Left panel: The differential production probability as a function of the azimuthal angle (j)e 
of the electron emission. Right panel: The anisotropy (j26p for the hyperbolic secant (solid curves) 
and Gaussian (dashed curves) shapes. For = 0.1 and C = 4. 

The differential probability of the e^e~ pair emission and the anisotropy as a function of 
the azimuthal angle 0e are exhibited in Fig. [51 The calculations are for the fast-decreasing 
one-parameter envelope functions for A = O.Svr, C = 4 and = 0.1. One can see a rapidly 
decreasing probability with 0e which leads to the strong anisotropy of electron (positron) 
emission. 

In case of the symmetrized Fermi distribution with small &/A, the situation changes 
drastically. As 6/ A —t- the envelope function goes to the fiat-top (step-like) shape /fs(0) 
^(A^ — (fP') with 9{x) = 1,0 for x > or x < 0, respectively, and correspondingly 

A 

y{l) = — [ ci0e*[''^-"''''('^-*''^] (27) 
27r J 

-A 

with / = / + $,'^(u. The function y{l) in the region ( < I < /max ^ 1 is alternating, rapidly 
oscillating with an amplitude that depends only on ^, (, and u. It is not sensitive to 0o- 
Modifications of 0o lead some phase shift of y{l) in a range of integration, leaving (|3^(0P) 
to be independent of 0o Therefore, the dependence of the integral of the partial probability 
w{l) ~ |3^(0P Eq. ( fT3l) on 0o is negligible. As an example, in the left panel of Fig. [6] 
we present the partial probability w{l) as a function of /, calculated at = 0.1, C = 4 and 
u = 1 for the small values of b/A equal to 0.1 and 0.01 at 00 = and vr, shown by solid 
and dashed curves, respectively. One can see some small modification of the frequency of 
oscillations at / ~ /min = C two extreme values of (po, but the amplitudes of the oscillations 
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FIG. 6: The partial probability w{l) defined in (jl3p at (pQ = and vr shown by solid and dashed 
curves, respectively, for the symmetrized Fermi envelope shape. The left panels correspond to 
small values of b/A =0.01 and 0.1, while the right panel is for b/ A = 0.5. For = 0.1 and C = 4. 

are similar. This situation is quite different from the case of the large value of 6/A = 0.5 
presented in the right panel of Fig. [6l One can see a strong difference in the /-dependence of 
w{l) for 00 = and tt. In the first case, the function w{l) has only one oscillation in a wide 
range of / and decreases smoothly with /. In the second case, the probability has a number 
of oscillations decreasing rapidly with increasing /. As a result, the total probability in the 
second case is much smaller. 

This behavior can also be understood from a different point of view. The integral over / 
of the derivative of the partial probability w{l) in Eq. ( IT3l) 

dw{l) d 



■\y{i\'-mi±i)y{i)\ 



(28) 



is vanishing because of the alternating and oscillating nature of y{l). Therefore, the proba- 
bility W is independent of 00 . 
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FIG. 7: The same as in Fig. [5] but for the symmetrized Fermi shape. The solid, dot-dashed and 
dashed curves are for b/A = 0.1, 0.3 and 0.5, respectively. For = 0.1 and C = 4. 
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In Fig. [7] we present our results for the symmetrized Fermi shape for the production 
probabihty (left panel) and for the anisotropy (right panel) for b/A = 0.1, 0.3 and 0.5. 
The result for b/A = 0.5 is similar to that shown in Fig. El However, for smaller values of 
b/A, the probability is a smooth function of 0o with some modest enhancement around 7r/2, 
which leads to a negligible anisotropy. 



IV. SHORT PULSES 



In this section we consider short pulses with the number of oscillation N > 2, however, 
many results are valid even for pulses with ~ 1. As we have seen, the one-parameter 
envelope shapes lead to similar results even for ultra-short pulses, therefore, later on we will 
limit our consideration to two extreme envelope shapes: hyperbolic secant and symmetrized 
Fermi shape with b/A = 0.1. 

As mentioned above, Eqs. f ll3p and (IT^ can be used for numerical estimates of the e~^e~ 
production probability evaluating five dimensional integral(s) with rapidly oscillating func- 
tions. Technically, such an approach needs long calculation time which makes it difficult for 
applications in transport/Monte Carlo codes. However, a closer inspection of the functions 
V{(t>) and C^*^(/) shows that the number of integrations may be reduced and, in some cases, 
Eq. ( fT4l) may be expressed in an analytical form. Thus, integrating by parts the function 
P(0) might be rewritten in the following form 



V{4>)^Vo{<P)-eCu I d<p'f\^'), Vo{<P) = z (^sin(0-0o)/(0) + o(^)) (29) 

— oo 

with 



O = "1 / - '^o)/(0') • (30) 

— oo 

The contribution of this term to V{(f)) is sub leading for the finite pulse size A = nN with 
N >2. First, because of the explicit factor 1/A, and second because the derivative /'(0) in 
the integrand reaches its maximum value at the boundaries of the pulse, where this function 
is suppressed. For an illustration, in Fig. |H] we present results of a numerical analysis of 
Vq{4') with the hyperbolic secant envelope function. The solid and dashed curves exhibit 
calculations with and without the term ( 130|1 . respectively for (pQ = and n. The left and 
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right panels correspond to A = ttA^ with N = 2 and 5, respectively. The term (|0(1/A)|) 
is shown by dot-dashed curves. One can see, in fact, that this term is rather small and may 
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FIG. 8: The function 'Po(<^) defined in (I29p with (soHd curves) and without (dashed curves) the 
term ()30p for A = vr with N = 2 and 5, shown in left and right panels, respectively. The 
term ()30p is shown separately by dot-dashed curves. 

be omitted. For the fiat-top envelopes this approximation is even much better. 

Using this approximation one can express the basic functions C^*^(/) defined in Eqs. (E]) 
and (fTT]) through the new functions Yi and Xi, which may be considered as an analog of the 
Bessel functions in IPA, 



-oo 
oo 



Xi{z) = — I #/2(^ + 0o)e^"^-^^^''^'^. 



r(0) 



f 



(31) 



where the function r{(j)) is a smooth function of (j). For the hyperbolic secant we have 
r{(f)) = Atanh(</)/A), where we skip the constant term which does not contribute; for the 
fiat-top envelope, r(0) ~ ^(A^ — 0^). The new representation of the basic functions C-*-* (/) 
reads 



1 



- (l^+ie*(^+^)<^° - l^.ie^C-i)-^") 
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= ll(^)e*Wo, Y,{z) = ^iY,^,{z)+Y,_,{z))-e-Xi{z) . (32) 

21 ui 

It allows to express w{l) in Eq. (HM in the form 

wil) = 2Y,\z) + e{2u - 1) (yU^) + r,2+i(^) - 2Yi{z)Xt{z)) , (33) 

which resembles the expression for the probability in case of IPA (cf. Eq. (lAip ). Now we 
are going to discuss separately the weak-, intermediate- and strong- field regimes. 

A. Production probability at small field intensities (^^ ^ 1) 

In case of small values of .^^ <^ 1, implying 2; < 1, we decompose I = n + where n 
is the integer part of /, yielding 

00 

Yi ^ ^ ! d^^'^-''''^^f{'4^ + (t>o) 
2tx J 

—00 

00 

= ^ / rf^5^^sin-^e^("+^)'^r+^(^ + 0o) . (34) 
271 J ^-^ ml 



-00 



m=0 



Similarly, for the function Xi{z) the substitution jm+2 g^ppiigg_ "Yhe dominant 

contribution to the integral with rapidly oscillating integrand comes from the term with 
m = n, which result in 



^n+e - ^e-<^«F('^+i)(e) , X„H., ^ ^e-^«F("+2)(,) , (35) 

where the function F^'^\e) is the Fourier transform of the function /^{ip). 

As an example, let us analyze the e~^e~ production near the threshold, i.e. C ~ 1- 
this case, the contribution with = 1 is dominant and, therefore, the functions Yo+e are 
crucial, including the first term in ( 133|) . The functions Xo+e are not important because they 
are multiplied by the small ^'^ and may be omitted. Negative e = ^ — 1 and positive e 
correspond to the above- and sub-threshold pair production, respectively. The function Fo+e 
reads lo+e = F^^\e) exp[— z0oe] with corresponding Fourier transform F^^\e) presented in 
Eq. ( 1221) . Note that the 0o-dependence of the production probability disappears in this case 
because the latter one is determined by the quadratic terms of the F-function. 
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Consider first the pair production above the threshold. Keeping the terms with leading 
power of C,"^ one can express the production probability as 



du 



47Vn 







Ml 





e 



u 



3/2^ 



(36) 



u 



where, taking into account that, at finite values of A, Fourier transforms for all considered 
envelopes decrease rapidly with increasing e one can get 

1/2 oo oo 

j deFW\e)c^ I rfeF«2(e) = i- J d<i> f\<P) ^ . (37) 

1~C — oo — oo 

Combining these two equations one recovers exactly the IPA result jsl. Thus, we can con- 
clude that for small field intensities for a finite pulse duration the probabilities of e^e~ pair 
emission above threshold with C < 1 results in a coincidence of IPA and FPA, independently 
of the shape of the envelope function. For an illustration, in Fig. [9] we show the partial 





FIG. 9: The partial probability w{l) defined in (133p as a function of / at u = 1. The solid and 
dashed curves correspond to the beam size A = vrA^ with N = 2 and 10, respectively. Left and 
right panels exhibit results for the envelopes with hyperbolic secant and symmetrized Fermi shapes, 
respectively. For = 10~^ and ( = 0.5 



probability w{l), calculated at n = 1 for the above-threshold region with C,"^ = 10~^ and 
( = 0.5 in a finite region of / for the envelope size A = vrA^ with N = 2 and 10, respectively. 
For the envelope with a hyperbolic secant shape (left panel) one can see smooth curves with 
maxima at integer values of /. The widths of bumps decrease with increasing N. However, 
the integral of w{l) over / in the neighborhood of the first maximum is independent of TV 
and coincides with the contribution of the first harmonic in IPA which leads to an equality 
of IPA and FPA results. For the symmetrized Fermi shape (right panel) the situation is 
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different in some sense. The corresponding Fourier transforms F^p (l) in (135|) oscillate with 

/. For example, the function F^p goes to zero at a multiple of This results in an 

oscillating structure of w{l). However, the exponential decrease of w{l) with increasing of 

the integer values of / is the same. 

The situation changes when we are slightly below threshold, i.e. C > 1. In this case, 

the function lo+e dominates again and the result for FPA is the same as in ( l36l) but with 

1 

the substitution Jo — /i, with h - J deF^^^^{e). In case of smooth envelope shape 

C-i 

(e.g. hyperbolic secant) the dominating contribution to this integral comes from the lower 
limit and, therefore, Ii ~ (^-^hPlC ~ l)j • As a result, the production probability strongly 
depends on the duration A of the pulse. 

In case of a fiat-top envelope, we have a similar effect, because Flp{l) in general decreases 
exponentially as exp(— vrW), where b increases with increasing at fixed b/A. 




FIG. 10: The same as in Fig. [9] but for the sub-threshold region at C = 1.25. 



In Fig. [iniwe show the partial probability w{l) in the sub-threshold region, i.e. ( = 1.25. 
One can see that for the hyperbolic secant envelope (left panel) the difference of w(/) at / ^ C 
for N = 2 and = 10 is more than several orders of magnitude, which will be reflected in 
the total probability. In the case of the symmetrized Fermi envelope shape, one also can see 
a significant enhancement of w{l) for N = 2 compared to = 10. But now, the difference 
between FPA and IPA is larger compared to the case of the hyperbolic secant shape. 

The total probability W of e^e~ emission as a function of the sub-threshold parameter ( in 
the vicinity C ~ 1 is presented in Fig. [TTl The left and right panels correspond the hyperbolic 
secant and symmetrized Fermi envelope shapes, respectively. Calculations are performed for 
short pulses with A^ = 2, 5 and 10 oscillations in the pulse at C,"^ = 10~^. For comparison, 
we present also the IPA results. In the above-threshold region, results of IPA and FPA are 
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FIG. 11: The total probability W of the e~^e~ pair production as a function of C for short pulses 
with A = vrA^ for = 2, 5 and 10 as indicated in the legend. The thin solid curves marked by dots 
depict the IPA result. Left and right panels correspond to the hyperbolic secant and symmetrized 
Fermi envelope shapes, respectively. 



equal to each other according to Eqs. fl36|) and fl37j) . However, in the sub-threshold region, 
where ( is close to unity, the probability of FPA considerably (by more than two orders of 
magnitude) exceeds the corresponding IPA result. In case of the hyperbolic secant envelope 
function the probability increases with decreasing pulse duration. The results of FPA and 
IPA become comparable at > 10. Qualitatively, this behavior is true for the case of the 
symmetrized Fermi distribution. However, in this case, the enhancement of the probability 
in FPA is much greater. This is due to the fact that the envelope of the maxima in the 
partial probability w{l) (cf. Fig. [TOl) decreases with increasing / in different ways for different 
envelope shapes. In case of the hyperbolic secant it decreases as exp(— vrA/), whereas in case 
of symmetrized Fermi shape it decreases as exp(— 27r6/). For the latter one, aX b/A = 0.1 the 
slope is much smaller. Such a strong gain of e~^e~ emission is expected for other values of ( 
when ( exceeds an integer number. This effect is illustrated in Fig. [121 where the total e'^e~ 
production probability W is presented in a wide region of C at ,^ = 0.01. For convenience, we 
show also results for two different pulse shapes simultaneously. For two oscillations in a pulse 
(left panel N = 2), for the hyperbolic secant shape one can see a regular enhancement of the 
probability W when ( exceeds the corresponding integer value. As a result, W{() in FPA 
is a smooth function, while a step-like dependence of the probability appears in IPA. For 
the flat-top, symmetrized Fermi distribution at C > 1, the probability is significantly larger 
than for hyperbolic secant pulse shape and displays a step-like behavior. The latter one, 
however, is related mainly to the oscillating nature of the corresponding Fourier transform 



21 




FIG. 12: The total probability W of the e^e~ pair production as a function of C for one- and 
two- parameter envelope shapes (dashed and solid curves are for hyperbolic secant and symmetrized 
Fermi shapes, respectively). The thin solid curves marked by dots depict the IPA result. Left and 
right panels correspond to the number of oscillation in a pulse N = 2 and 10, respectively. 

in dHSD. 

At large values of (right panel, = 10) results of FPA and IPA become close to each 
other, especially for the one-parameter envelope shapes. For this case, at least for C, = 0.01, 
A^ ~ 10 can be considered to be near infinite, when considering the overall ( dependence. 
For the flat-top shape with small b/A the probability in FPA is higher than the result of 
IPA near integer values of (. 

To summaries this part we have to note that temporal beam shape effects for short pulses 
are strong and even dominant at small field intensities in the parameter region where the 
variable z is small, z <^ 1. At finite z, the non- linear dynamics of in the strong pulse 
becomes essential. 

B. Production probability at intermediate field intensities (^^ ~ 1) 

At finite values of z, z ^ 1, the probability of e~^e~ emission needs to be calculated 
numerically using Eqs. (fT3|) . ( 13T]) and (!33|) . In Fig. [131 we present the total probability W 
as a function of ( at fixed = 1 (left panel) and as a function of at fixed C = 4 (right 
panel). The calculations are performed for the hyperbolic secant and symmetrized Fermi 
pulse envelope shapes, shown by the dashed and solid curves, respectively. The duration 
of the pulse is A = vrA^ with N = 2. For comparison, we also present IPA results by the 
thin solid curves marked by dots. At finite the probability decreases monotonically with 
increasing ( (left panel), contrary to the step-like decrease typical for the small .^^ <^ 1 (cf. 
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FIG. 13: The total probability of e~^e~-pair production for two envelope shapes (dashed and 
solid curves are for hyperbolic secant and symmetrized Fermi shapes, respectively). The thin solid 
curves marked by dots are the result of IPA. Left panel: The total probability as a function of ^ 
at = 1. Right panel: The total probability as a function of at = 4. 

Fig. [T2l) . The probability for the flat-top pulse shape slightly exceeds the probability for the 
hyperbolic secant and the IPA result. 

Concerning the dependence (right panel), one can see a significant enhancement of the 
total probability W at small values of for the fiat-top pulse shape compared to the case 
of hyperbolic secant and the IPA result. The latter two results are practically identical to 
each other. At .^^ > 1, the production probability does not sensitively depend on the pulse 
shape, and FPA and IPA results are close to each other. This means that at large field 
intensity the dynamical aspects of the pair production gain a dominant role in comparison 
with the pulse shape and size effects. 

Finally, we note that, at finite the dependence of the probability on the azimuthal 
angle 0e disappears and the distribution in the x — y plane becomes isotropic. As an 

10^ 

I 

^ 10^ 

i.o- 

(N 10 

10^0 45 90 135 180 

(l>e(deg) 

FIG. 14: The differential probability of e'''e~-pair production as a function of 0e = (/"o at (" = 4 
and N = 2. 
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example, in Fig. [T3]we present results of calculations of the differential probability of e'^e~- 
pair production as a function of 0e = 0o at C = 4 for the hyperbolic secant pulse shape 
with = 2 at = 0.1, 1 and 10. The results reflect the isotropy of the e~^e~ emission and 
expose the dependence. 



C. Production probability at large field intensity (^^ » 1) 

At large values of 1, the basic functions Yj and Xi in Eq. (PTI) can be expressed 

in the form of flTH]) : 

oo oo 

Yi= J dqF^'\q)G{l-q) , Xi = j dq F^'\q) G{1 - q) , (38) 

— oo — oo 

where F^^\q) and F^'^\q) are Fourier transforms of the functions /(</>) and (</>), respec- 
tively, and G{1) may be written as 

oo 

^(0 = — / c/^e^'^-^'^'^^+^'M . (39) 
2tx J 

— oo 

In deriving this equation we have considered the following facts: (i) at large the probability 
is isotropic, therefore we put 0o = 0, (ii) the dominant contribution to the rapidly oscillating 
exponent comes from the region 0^0, where the difference of two large values 1(f) and z sin 
is minimal, and therefore, one can decompose the last term in the function V{(j)) in ( 129|1 
around = 0, and (iii) replace in exponent /(0) by /(O) = 1. 

Equation fl39l) represent an asymptotic form of the Bessel functions 3J] Ji{z) with / = 
I + S,'^Qu at 1, z ^ 1, and therefore the following identities are valid 

G'(/-1)-G'(/ + 1) = 2G;(/), G'(/-l) + G(/ + l) = 2-G(/) , (40) 

z 

which allow to express the partial probability w{l) in (l33l) as a sum of the diagonal (relative 
to /) terms: Yj^, YiXi, Xf and 1^'^. The integral over / from the diagonal term can be 
expressed as 

oo oo 

Iyy = j dlY^ = j dqdq'F^^\q)F^^\q') j dlG{l - q)G{l - q') . (41) 
C C 
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Taking into account that for the rapidly oscillating G functions G{1 — q)G{l — q') ~ S{q 
q')G^{l - q) and (g) < (/) ~ one gets 



1 

2^ 



^YY I dlG\l) 
C 



(42) 



d(t)f{(t)) J dlG\l) = Ny 

— oo ^ 

Similar expressions are valid for the other diagonal terms with own normalization factors. 

oo oo 

For the Xf term it is Nxx = ^ J d(Pf\^), and for Y^X,, Nyx = ^ J d^f\(p). At large 

— oo — oo 

the probability does not depend on the envelope shape, because only the central part 
of the envelope is important. Therefore, for simplicity, we choose the fiat-top shape with 
Nyy = Nyx = Nxx = Nq = A/it which is valid for any smooth (at ~ 0) envelopes. 

Making a change of the variable 1-^1 = 1 + ^"^Qu the variable z takes the following form 



4e'C' {uui 



UUr 



(43) 



with Iq = ({1 + and uj = 1/Iq, that is exactly the same as the variable z in IPA with 
the substitution I ^ I. All these transformations allow to express the total probability in a 
form similar to the probability in IPA for large values of ^'^ and a large number of partial 
harmonics n, replacing sum of n by an integral over n {g] 



oo "-i 

W = ^aM,C/^ j dlj 



du 



lo 



jf{z)+ei2u-i] 



' z^ 



{--l)jf{z) + J%z) 



.(44) 



Utilizing Watson's representation 34| for the Bessel functions at /, 2 ^ 1 and / > z, 
Ji{z) = (27r/ tanh a)"-'^/^ exp[— /(a — tanh a)] with cosh a = l/z, and employing a saddle point 
approximation in the integration in (144|) we find the total probability of e"'"e" production as 
(for details see Appendix A) 



3 /3aMee , 
^=8V2TV^^"^P 



15^ 



6C V 8C 



(45) 



This expression coincides with the production probability in IPA which is the consequence 
of the fact that, at ^ 1 in a short pulse, only the central part of the envelope at ^ 
is important. Approximating d = 1 + 0{^/Q, the leading order term recovers the Ritus 



result 



6|- 



For completeness, in Fig. [15] (left panel) we present FPA results for not too large values 
of C,"^ calculated for the hyperbolic secant envelope shape with N = 2 (curves are marked 
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FIG. 15: The total probability W of the e~^e~ pair production as a function of for various 
values of C- Left panel: Results of FPA for not too large values of (curves marked by "stars" 
in "FPA" sections) and the asymptotic probability ()45p for large values of (sections labeled by 
"asymptotic") at = 2, 4 and 6. Right panel: The asymptotic probability (jl5]) for various values 
of C as indicated in the legend. 



by "stars") and the asymptotic probabiUty calculated by Eq. ( l45l) at C = 2, 4 and 6. The 
transition region between the two regimes is in the neighborhood of ^ 10. In the right 
panel, we show the production probability at asymptotically large values of for 5 < C < 20. 
The exponential factor in psl) is most important at the relatively low values of ^"^ ~ 10 (large 
C/O- extremely large values of (small C/C ) the per-exponential factor is dominant. 



V. SUMMARY 



In summary we have considered different aspects of e~^e~ pair production in a strong 
electromagnetic field of a finite (laser) pulse, thus generalizing the Breit- Wheeler process to 
non-linear (i.e. multi-photon) effects. The pair production in the sub-threshold region with 
C > 1 is currently a subject of great interest. We have shown that the production probability 
is determined by a non-trivial interplay of two dynamic effects. The first one is related to 
the shape and duration of the pulse. The second one is the non-linear dynamics of charged 
particles in the strong electromagnetic field itself, independently of the pulse geometry. 

These two effects play quite different roles in two limiting cases. 

(i) The pulse shape effects are manifest clearly at small values of product the where ^ 
characterizes the laser intensity and ( refers to the threshold kinematics. The rapid variation 
of the e.m. field in a very short pulse amplifies the multi-photon events, and moreover. 
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the probability of multi-photon events in FPA can exceed the IPA prediction by orders of 
magnitude. Thus, for example in case of an ultra-short (sub cycle) pulse with the number 
of oscillations in the pulse less than one, the production probability as a function of ( is 
almost completely determined by the square of the Fourier transform of the pulse envelope 
function. High-/ components, where / is the Fourier conjugate to the invariant phase variable 
0, lead to the enhancement of the production probability. Among the considered envelope 
shapes, the flat-top shape with small 6/ A is most promising to obtain the highest probability. 
We also find that the different envelope shapes lead to anisotropics of the electron (positron) 
emission which can be studied experimentally. For short pulses with A^ < 10, the effects 
of the pulse shape are also important and the final yield differs significantly from the IPA 
prediction. This difference depends on the envelope shapes and the pulse duration. 

(ii) Contrary to that, the non-linear multi-photon dynamics of in the strong electro- 
magnetic field plays the determining role at large field intensities, ^ 1. Here, the effects 
of the pulse shape and duration disappear since the dominant contribution comes from the 
central part of the envelope function. As a result, the probabilities in FPA and IPA coincide. 

In the transition region of intermediate intensities ~ 1, the probability is determined 
by the interplay of the both effects, and they must be taken into account simultaneously 
by a direct numerical evaluation of the multi-dimensional integrals with rapidly oscillating 
integrands. 

Finally, we emphasize that the elaborated methods can be applied easily in transport 
approaches aimed at studying e^e~ pair production in the interaction of electrons/positrons 
and/or photons with a finite electromagnetic (laser) pulse. 



Acknowledgment s 

The authors acknowledge fruitful discussions with H. Ruhr, D. Seipt, and T. Nousch. 
The support by R. Sauerbrey and T. E. Cowan is gratefully acknowledged. 



Appendix A: Production probability at large values of ^ 

The total probability W in the limit of large ^ and and small C,/C, was evaluated by 



Narozhny, Nikishov and Ritus jo]. for completeness and easy reference, we recall here some 
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details of evaluation making expansion for an arbitrary ^/(. 

In IPA, the total probability is represented as an infinite sum of partial harmonics 



1 r rl 

-aM^C J ^,/,J^ _ ^ {^Jl{z) + i\^n - 1) {Jl^,{z) + Jl_,{z) - 2Jl{z)) } , 

(Al) 



^ 4 

n=no ^ 



where rio = rimin = ({1 + C,^), Un = n/uQ, and Jn{z) is the Bessel function of the first kind 
(cylindrical harmonics). Using the identities 

Tl 

2 - Jn{z) = Jn-l{z) + J„+l(2), 2 j' n{z) = Jn-l{z) - Jn+l{z) , (A2) 
Z 

the total probability takes the following form 
W = \aM^C^^Y. I ^3/2%^ (^n(^) + ei2u - 1) (^(g - imz) + J'liz)^ ^ .(A3) 

At large ^^1,(^^1,72, 2;^! and n > z one can replace the sum over integer n by an 
integral over dn, replacing, for convenience, integer n to continues I with /mm = = C(1+'C^)- 
Using Watson's asymptotic expression for the Bessel functions one finds 

Jl ( = , ^ g-/(»-tanha) ^ Q (\\ (A4) 

V cosh ay V27r/tanha \i) 
with cosh a = 11 z. If / is large the first term represents a good approximation irrespectively 



whether ^/C is small or large [3J| . The corresponding derivative reads 

J[[z) ~ sinha Ji(z) (1 + ^ ) . (A5) 

\ 2/ sinh a tanh a J 

Consider first the case of small ^/C <^ 1, when the second term in (lASP can be neglected. 

Then, the total probability becomes 

00 Ul 

e'M^C'^ f f du l + 2e^(2n-l)sinh^a 

W = — — — / dl / — — r exp / u,/ , A6 

o-w^ J J u^/^yju — l / tanha 

lo 1 

where ui = 1/Iq and f{u, I) = —2l{a — tanh(a)) with 

1+^2 (i_2uy 

tanh^(a) = ^ ^ . (A7) 

To avoid a notational confusion with respect to the standard variable a, we replace below 
the fine structure constant by e'^/in. 
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The two-dimensional integral is evaluated using the saddle point approximation since 
the function f{u,l) has a sharp minimum at the point u = u defined by the equation 
f^{u = u) = 0. That allows (i) to expand it to a Taylor series 

f{uj) ^ + n)2 , (A8) 

and (ii) to take the rest (smooth) part of the integrand in Eq. ( 1A6I) at the point u = u 
yielding 

OO Ul 

167r^ J J Vu-l 



with 



The explicit expression 



leads to the solution 



2 2/o 

which results in the following equalities 



p.. 4/osinh a ^'^ / 2u\ 



S = ^ = ;^, (A12) 



2 1 -~ 8/ 

tanha = tanha(n) = — ^ sinha = -, 



A = ^^^^Sr^A/lTe^, /(«,/) = -2Z(a-tanha) . (A13) 

Using the substitutions m = t + 1, a = 2(a — tanha), and A = —^f"(u, I) one can rewrite 
Eq. ^M) as 



OO OO 

/.M,(M)e-— 7.*r-.e-'-. ,AM, 



167r2 

lo 

with z/ = 1/2, /3 = y4, and 7 = 2^4(1 — u). The integral over dt is expressed via the parabolic 
cylinder function D^^ 

* <-'e-^'-- = (i^) r(.) explgl D,. (-2^) . (A15) 
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which results in 



167r3/2 



■ / dl ' A(n,/)e-'^'"4(i--)^Z}„.(y) (A16) 



with y = v2A(l — u). The main contribution to this integral comes from the region m ~ 1 
(/ ~ r = 2/0) and, therefore, one can use the substitution 



00 



dl = -^L [ dy^^ I dy, (A17) 



2A J V2A 

'0 Ja/2 



which results in 



W 

167r3/2 

with Z = 2loa/V2A. Using the identity 



00 

/ 



d|/e^^-^'/^D_i(l/) = J^e^'/^ ^ (A19) 
2 \ Z 



one can rewrite the production probability as 

W='-^^^^Ao{uJ)exp[-2loa+'-^] . (A20) 

In order to reproduce the Ritus result 6| in terms of the kinematic factor C and the field 

intensity C, one has to use the identity /q = C(l + ■C^) ^'^^ to represent a{a) as a series for 

small values 1/^ utilizing the expansions 

,111 3 , _ 1 11 3 , 3 

a = arsmh- ~ — + —— ^ , tanha = — , ~ — H —, Ar, = —— A21) 

^ ^ 6^3 40^5 ' ^/TT^ ^ 2^3 8^5' 'J 2C^ ^ 

which leads to ( l45l) with (i = 1. Inclusion of the second term in ( lASp modifies eventually ^0 

as 

^ = ^7^1 + 7^^ + 1:11 (A22) 



2Ce V 6C V 8C. 

yielding the result displayed in ( 145|) which extends the Ritus result for arbitrary values of 
^/(. We emphasize that, in the strong field regime, IPA is representative since, as stressed 
above, pulse shape and pulse duration effects are sub leading. 
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